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Abstract 

We compute the G(a^) two-loop corrections to the neutral CP-even Higgs boson mass ma¬ 
trix in the Minimal Supersymmetric Standard Model, for arbitrary values of m A and of the 
parameters in the stop sector, in the effective potential approach. In a large region of pa¬ 
rameter space these corrections are sizeable, increasing the prediction for nih by several GeV. 
We present explicit analytical formulae for a simplified case. We discuss the inclusion of 
momentum-dependent corrections and some possible ways of assigning the input parameters. 


1 Introduction 


The Minimal Supersymmetric extension of the Standard Model, or MSSM (for a recent review 
and references, see e.g. Ref. [|1J]) is today a paradigm for possible new physics at the weak scale, 
supported by many theoretical and phenomenological considerations. The predictive power of 
the MSSM, however, is strongly limited by our ignorance of the mechanism for supersymmetry 
breaking in the underlying fundamental theory. In particular, the spectrum of supersymmetric 
particles, with its strong dependence on many soft supersymmetry-breaking parameters, is es¬ 
sentially unconstrained. This is not true for the MSSM Higgs sector: the dimensionless Higgs 
couplings in the scalar potential are controlled by supersymmetric Ward identities that receive 
only finite quantum corrections from the soft breaking of supersymmetry. As a result, the masses 
and the couplings of the MSSM Higgs bosons, (h, H, A, H^), can be expressed, at the classical 
level, in terms of the gauge boson masses (m w ,m z ) and only two additional parameters: for 
example, the mass of the neutral CP-odd state, m A , and the ratio of the two Higgs vacuum ex¬ 
pectation values (VEVs), tan (3. The number of independent predictions in the Higgs sector does 
survive at the quantum level, but radiative corrections induce finite modifications of such predic¬ 
tions, and bring along additional dependences on the remaining MSSM parameters, associated 
with the virtual particles in the loops. This is the reason why an impressive theoretical effort 
has been devoted to the computation of the radiative corrections to the MSSM Higgs sector. 

In this paper we address once more the crucial issue of the computation of the neutral CP 
even Higgs boson masses. Soon after the first calculations of the dominant one-loop effects, 
associated with top and stop loops [|j], the full one-loop computation became available |jj, f|, 
and resummations of leading and next-to-leading logarithms were included [|6| via appropriate 
renormalization group equations (RGEs). We are now at the stage in which the computation 
of the most important two-loop effects is being completed. In analogy with the one-loop case, 
the most important two-loop effects are generated by strong and top-Yukawa corrections to 
the one-loop diagrams involving the top quark and its supersymmetric partners, i.e. they are 
0(ata s ) and 0(a 2 ), respectively. 

The first genuine two-loop computation [pj considered 0(ata s ) and 0(a 2 ) corrections to m^, 
and was restricted to the limiting case of tan (3 —* oo and negligible stop mixing. 

At 0(ata s ), the result of @ was generalized by three groups, all working in the limit of vanish¬ 
ing external momentum, i.e. in the effective potential approximation (EPA). The diagrammatic 
computation of [|| is valid for arbitrary values of m A and of the stop parameters, but in such 
general case its complete formulae are rather lengthy and available only in the form of a com¬ 
puter code. The computation of || is applicable to the corrections to m^, in the limit m A m z , 
and its results are available in analytic form under further restrictions. The computation of [|1C]] 
has recently provided simple analytical formulae for the two-loop corrected mass matrix of the 
neutral CP-even Higgs sector, for arbitrary values of m A and of the stop parameters. 

The inclusion of the 0(a 2 ) corrections, which can be of comparable numerical importance, 
is still at an earlier stage. After |t|], the only genuine two-loop computation so far is the one 
of [11], which, as the one of Q, is applicable to the corrections to rrih, in the limit m A m z . 
Moreover, the full results of [pd]] are available only in numerical form, but simple analytical 
formulae are provided under the additional assumptions of large and universal soft stop masses, 
m 2 Q = mfj m 2 . In the present paper we go one step further, and compute the O(af) two-loop 
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corrections to the mass matrix of the neutral CP-even Higgs bosons. We still work in the EPA, 
but our calculation is valid for arbitrary values of m A and of the stop parameters. 

Our paper is organized as follows. After this introduction, Section 2 recalls some general 
features of the calculation of the neutral Higgs boson masses in the MSSM. In particular, we 
describe the connection of the EPA with the full calculation, and an ‘improved’ perturbative 
determination of the masses from the full momentum-dependent two-point function. Section 3 
describes the main features of our two-loop calculation of the 0(a 1 2 ) contributions. Explicit 
analytical formulae are displayed for the simple case of large and universal soft stop masses, 
but arbitrary values of all the remaining Higgs and stop parameters. Section 4 discusses in 
detail some possible ways of assigning the input parameters, and the relation between the DR 
scheme and on-shell (OS) renormalization schemes, identifying a version of the latter that we find 
particularly convenient. Section 5 presents some quantitative evaluation of the 0(a 2 ) corrections 
to rrih and m H , and numerical comparisons with the previously available results, for a number 
of representative cases. We find that the 0(a 2 ) corrections can be sizeable in a large region of 
parameter space, increasing the prediction for m/,, by several GeV. In particular, for large stop 
mixing most of these corrections are genuine two-loop effects, which cannot be accounted for 
by standard renormalization-group improvements. In the concluding section we summarize our 
results and comment on possible extensions. Technical details, such as the explicit formulae that 
are needed for the transition from the DR scheme to our implementation of the OS scheme, 
are confined to Appendix A. We find that our general analytical result is too long even for an 
appendix, thus we make it available, upon request, in the form of a computer code Q 


2 General results 


We review here, along the lines of Q, some general results that apply to the calculation of the 
MSSM neutral Higgs boson masses, in the CP-conserving case, at every order in the different 
coupling constants and in the loop expansion. 

In the EPA, the CP-odd and CP-even mass matrices are identified with the second derivatives 
of the effective potential, 14s = Vq + V, evaluated at its minimum: 


/ 9 \eff 

44 - 


d 2 V P 


eff 


8P i dP j 


/ 9 \ eff 

44 - 


d 2 14 


eff 


ij dSidSj 


(i,j = 1,2), 


(1) 


where we have decomposed the Higgs fields into their VEVs plus their CP-even and CP-odd 
fluctuations as follows: 


0 _ vi + Si + iPi 






V2 + 5*2 + i P r 2 


( 2 ) 


42 ’ v/2 

It is understood that I4ff is expressed in terms of renormalized fields and parameters, in the 
Landau gauge and in the DR renormalization scheme In particular: 

k eff / n \ 0, eff / n \ eff 


(Mi) - 4l) + (AVI) . 


( 3 ) 


1 E-mail: pietro . slavichOpd. infn. it. 

2 In our decomposition, Vo has the same functional form of the tree-level potential, whereas V contains the 
residual loop corrections. The gauge choice will not be relevant for the two-loop calculations of the present paper. 
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where 


( m l $ + s 2 p ~ (m\ + m 2 A ) sp cp \ 

V ~(m 2 z + ml) Sp Cp m 2 z s 2 p + m 2 A c 2 p J ’ 



— (—l) i+J 


d 2 V 

dPjdPj 


+ 

min 


d 2 V 

dSidSj 


min 


( 4 ) 

( 5 ) 


In Eq. (||), m 2 A is the mass of the CP-odd state in the EPA, i.e. the non-vanishing eigenvalue 
of (ASimilarly, = ( g 2 + g' 2 )v 2 /4 is just the DR mass for the Z boson Q, where 
v 2 = v 2 + v 2 , and v\ and v 2 are the DR VEVs computed by minimizing I4ff- Finally, sp = sin (3 
and cp = cos/3, where tan/3 = V 2 /V 1 . 

Since I4ff generates one-particle-irreducible (1PI) Green’s functions at vanishing external 
momentum, it is clear that the EPA neglects some momentum-dependent effects in the Higgs 
self-energies, and the entire finite contribution to the Z self-energy. The complete computation 
requires the full two-point function [] T s(p 2 )i which takes the form of a 2 X 2 matrix in the space 
(Si, S 2 ), and the analogous two-point functions for the neutral CP-odd spin-0 bosons and gauge 
bosons. Using the standard decomposition of the two-point function, 


r s(p 2 ) = P 2 - M 2 s {p 2 ), 


( 6 ) 


the physical masses (ml,m 2 H ) are the two solutions of the equation 


det 


P 2 -M 2 s (p 2 ) 


= 0. 


( 7 ) 


At any loop order, a full diagrammatic calculation relates these two masses with the physical 
masses of the neutral CP-odd Higgs boson ( m 2 A ) and of the Z vector boson (m|), defined in a 
similar way. 

To make contact between the complete calculation and the EPA, we need the relations 
between [M. 2 s (p 2 ), m 2 A1 m 2 \ on the one side, and [(Aon the other side. The 
transverse ZZ propagator and the A A propagator have the form [p 2 — m| + n^^(p 2 )] -1 and 
[p 2 — ml + respectively. Thus we have 


ml = m 2 z + Ii Z z(m 2 z ), (8) 

ml = m 2 A + 5U AA (m 2 A ). (9) 


It is understood here that the (DR-renormalized) two point functions n zz(p 2 ) and 5Haa(p 2 ) 
include both 1PI and non-lPI (mixing) terms. Indeed, beyond the tree level the transverse 
Z boson mixes with the photon, whereas the the field A = sin (3 P\ + cos /3 P 2 mixes with the 
Goldstone boson G = — cos (3 P\ + sin/3P2 and with the longitudinal Z boson. In the gaugeless 
limit, for instance, 5Haa(p 2 ) reads 

_ m ^ 2 >= An -<^ - JTaS) • < io > 

3 This choice is not imposed by the EPA, but is rather a matter of convenience. Alternatively, we could include 
in m| the self-energy corrections evaluated at vanishing external momentum. 

4 We absorb in our P’s the usual i factor and consider only their real part. The same conventions apply then to 
the self-energies. 
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( 11 ) 


( 12 ) 


where, for a generic self-energy IT a b(p 2 )i we define 

All ab (p 2 ) = U a b(p 2 ) - n afe (o). 

For the two-point function of the CP-even Higgs fields, we can write 

Ml ( p 2 ) = (Ml)° + (A M%f + (AMlf , 

where (Mg)° has the same form of (,A/f;|) 0 ’ eff in Eq. (|4j), with the only difference that the EPA 
masses (m A ,m |) are replaced by the physical masses (m A ,m 2 ). The matrix (AAl|) eff was given 
in Eq. (!)• Finally, the matrix (AMg ) p2 has entries 

n Z z(m 2 z ) c 2 0 + 5U AA (m 2 A ) s 2 0 - An 5lSl ( p 2 ), 

- n zz (m|) + m AA (m 2 A ) s 0 ci 3 - An 5l 5 2 ( p 2 ), 
n zz{m 2 z ) 4 + 5n AA (m 2 A ) 4 - An S252 ( p 2 ). 


AMl) 


( 

(a.m|) 


( 


AM 


12 
2\P 2 


S ) 


22 


(13) 

(14) 

(15) 


This establishes the desired correspondence between the EPA and the full computation. 

We now describe how to perform an ‘improved’ perturbative determination of the neutral 
CP-even Higgs masses (and of the associated mixing angle), starting from Eq. ©■ Of course, 
this equation could be solved numerically, at the cost of some computer time. Our solution 
has the advantage of producing explicit and accurate analytical formulae, which can facilitate 
a number of checks and considerably speed up the numerical evaluation of (m|, m 2 H ), for any 
given set of input parameters. The latter feature is welcome for the practical implementation of 
our results in computer codes. We have explicitly checked that, for a wide range of parameter 
choices, our perturbative solutions of Eq. 0 agree with numerical solutions within 10 3 GeV. 

The matrix M 2 - ( p 2 ) can be conventionally decomposed into the sum of a p 2 -independent 
‘unperturbed’ piece and a p 2 -dependent ‘perturbation’: 


M 2 s (p 2 ) = AT + AM 2 (p 2 ) 


(16) 


At the lowest order, i.e. neglecting AM 2 (p 2 ), Eq. ( 0 ) reduces to the diagonalization of M . This 
gives the eigenvalues 


m H,h 


1 

-9 -9 

1 ( - 9 -2 \ 2 

/-9 \ 2 

2 

+ A^22 \j 


+4(a*; 2 ) 


and the eigenstates 



cos a sm a 
— sin a cos a 



where the mixing angle a is determined by 

t 2 


sin 2 a = 


2 M 


12 


rn H — m 


2 > 


cos 2 a = 


M 2 n -M 2 


22 


mi, — mi 


(—7t/2 < a < 7r/2) 


At first order in AM 2 (p 2 ), we can write for the two eigenvalues 

[m 2 ] 1 = + sin 2 a AMh(mf l ) — sin2cf A_A4? 2 0^4) + cos 2 a AM\ 2 {jn\) 


[m 2 H ] 1 = m z H + cos^ a AMii(m' z H ) + sin 2 a AATf 2 ( r?1 iff) + sin^ a AM^i^H) 


—2 


2 — 


a 2 t —2 


rrr A A (—?■ 


,2 — 


<2 


(17) 

(18) 

(19) 

( 20 ) 
( 21 ) 
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At the second order: 


m h = m h + sin 2 a AAdi 1 ([m|] 1 ) — sin 2a AMf 2 ([mlY) + cos 2 a A -^22 


2il\ 


,2nl > 


_9 _9 

mir — mi 


sin 2a 


AMi 2 (m 2 h ) - AM' 2 n (m 2 h ) + cos 2a A M{ 2 (m 2 h ) 


( 22 ) 


m 2 H = m 2 H + cos" a AM.\ l {[m' H } ) + sin 2a AAdj 2 ([ m «] ) + sin 2 a AA4 22 ([m 2 H \ ) 


+ 


m\ - m 2 h 


sin 2a r 


AM' 2 2 (m 2 H ) — AAdfi (m 2 H ) + cos 2a AM.\ 2 (m 2 H ) 


(23) 


As discussed in [||, the decomposition of Eq. ([id]) is not unique. It is convenient to specify it 
in a way that improves, as much as possible, the accuracy of the perturbative solutions of Eq. ( 0 ), 
given by Eqs. (22) and (23), and allows to absorb an important part of the corrections to the 
Higgs interaction vertices into the definition of the renormalized mixing angle a. The natural 
choice M." = (Ad|)° and A M. 2 (p 2 ) = (AAf|) eff + (A M. 2 s ) p2 can work well for m A m z , 
while it is known H not to be very effective for m A ~ m z . A choice that leads to an improved 
convergence, and will be adopted in our numerical evaluation, is instead 

. 0 / . . n \ eff ' - - 2 


M 2 = ( Ml ) + (am|) C , A M 2 (p 2 ) = ( AM 2 s y 


(24) 


although other choices can be equally valid. 

Once the decomposition of M? s (p 2 ) has been specified, the perturbative solutions described 
above can be written more explicitly, taking also into account the conventional loop expansion 
of the self-energies (n = + ...). Consider the decomposition specified by ([h|). Then 

the two-loop corrected Higgs masses can be written as 

m h = , (25) 

m 2 H = m 2 H + S^rriff + 5^m 2 H , (26) 

where the two-loop corrected masses obtained in the EPA, (m 2 , fn 2 H ), are to be supplemented 
by the corrections that stem from (AM 2 s ) p2 . At first order we have (3|: 

S^ml = sin 2 ((3 + a) n {2 l(m 2 z ) +cos 2 (/3-a) AU { }\(m 2 A )-An^(m 2 h ), 


5 {1) m 2 H = cos 2 (/3 + a)n§^(m|)+ sin 2 (/3-a) An^(m 2 ) - An^(m|) 


(i) 




(27) 

(28) 


At the next order: 


5 {2) m 2 h 


= sin 2 (P + a) U {2) z (m 2 ) + cos 2 (P - a) <HI^(m 2 ) - AU^(m 2 h ) 


xtt( 2 ) 


r (2)/_2. 


Au hh( m l + 5[1)m l) - An S( m i) 


m 2 H - m 2 h L 


siii 2(P+a) ^(i) 


sin 2 (P~a) 


u zz( m l) + 0111 ^ An AA (™ 2 a)+ An S( ml) 


, (29) 


^( 2)^.2 


m„ = 


+ 


cos 2 (P + a) n^(m 2 ) + sin 2 (P - a) SU ( A ’ A (rn 2 A ) - AU { ^ H (m 2 H ) 
AU {2) H (m 2 H + 5 {1) m 2 H ) - An^(m|) 


1 


rn 2 H — rn'l 


sin 2 (P+a) ^(i) 
zz 


n (m 2 z ) + Sm2 ^f ft) AnJJ^ (m 2 A )+AU { ^ h (m 2 H ) 


1 2 


(30) 
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The corrections (S^rrif^, 8^m%) involve one-loop self-energies that are known and should be in¬ 
cluded to correct the EPA results ||]. On the other hand, the two-loop effects in <y( 2 ) m? H ) 

arise either from products (or iterations) of one-loop self-energies, or from 1PI two-loop con- 

(o) (o') (O') (0) 

tributions in (Jl K zz ,8Tl\ A , /^Ir HH ), which have not been computed so far. However, the 

effect of S^m 2 is always sub leading, and the same is true for 8^m% when itla is not large. 

To qualify the latter statement, we recall that the dominant two-loop corrections to (m\,m 2 H ) 
are 0(ata s m 2 ) and 0(a 2 m 2 ), described elsewhere as 0(ata s + a 2 ) for brevity. In particular, 
contributions of that order are present in the two-loop corrected EPA masses (m?, Wi 2 H ). Also 
, 8 ^m%) give 0(a 2 m 2 ) terms, induced by the one-loop corrected ( mf l ,m 2 H ,a ). The 
question is whether ffi 2 + S^rri^ and m 2 H + 8 ^m% exhaust the full 0 (a t a s m 2 ) and 0 (a 2 m 2 ) 
corrections to the Higgs masses or ( 8 ^m%, 8 ^m%) give additional contributions of the same 
order. When m 2 A = 0(m 2 ), the contributions from ( 8 ^m\, m 2 H ) are at most 0(atot s m 2 z ) or 
0(a 2 m 2 z ), i.e. they are suppressed by 0(m 2 /m 2 ) with respect to the dominant terms. In the case 
of 8 ^ 171 %, this is true also when m 2 A > m 2 z . Indeed, the terms growing with m 2 A in 8 H^%(jn A ) or 
An^(rr^) are actually suppressed by the associated coefficients, since cos 2 (j3 — a) = 0(m 4 z /rn 4 A ) 
and sin2(/3 — a) = 0(m 2 z /m 2 A ). Therefore the 0(a t a s m 2 ) and 0(a 2 m 2 ) corrections to m\ 
are fully accounted for by + S^m^: no further two-loop calculations are needed, as first 
emphasized in 0. On the other hand, 0(ata s m 2 ) and 0(a 2 m 2 ) contributions do generically 
appear in 8 ^m% when m 2 A > m 2 z . In the window m| < m 2 A < m 2 , however, such effects are 
suppressed at least by a factor 0 (rri A /rn(). 

3 0(af) two—loop corrections to the neutral Higgs boson masses 

We shall now describe our two-loop computation of the matrix (AAl|) eff at 0(a 2 ), for arbitrary 
values of the relevant MSSM input parameters. The computation is consistently performed 
by setting to zero all the gauge couplings (the so-called gaugeless limit), and by keeping h t = 
y/‘i-nott as the only non-vanishing Yukawa coupling. In this limit, the tree-level (field-dependent) 
spectrum of the MSSM simplifies considerably: gauginos and Higgsinos do not mix; charged and 
neutral Higgsinos combine into Dirac spinors with degenerate mass eigenvalues, |//| 2 , where /r 
is the Higgs mass parameter in the superpotential; the only massive SM fermion is the top 
quark; all other fermions and gauge bosons have vanishing masses; besides the stop squarks, 
to be discussed in detail below, the only sfermion with non-vanishing couplings is bi, a mass 
eigenstate with eigenvalue mi = mg; the lighter CP-even Higgs boson, h, is massless, as can 
be seen immediately from (||) with m| = 0; the same is true for the Goldstone bosons; all the 
remaining Higgs states, ( H , A, H^), have degenerate mass eigenvalues m\. The tree-level mixing 
angle in the CP-even sector is just a = f3 — ir/2. 

According to Eq. (||), (AM|) eff can be computed by taking the derivatives of V with respect 
to the CP-even and CP-odd fields, evaluated at the minimum of I4ff- In this computation 
we follow closely the strategy developed in [|10|], and express the field-dependent masses and 
interaction vertices that contribute to V, and are relevant for our calculation, in terms of five 
field-dependent quantities, which can be chosen as follows. The first three quantities are the 
squared masses of the top quark and squarks, 

m 2 = h 2 \H° 2 \ 2 , (31) 
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(32) 


mi = 
1 1,2 


In the above equation, 


(m 2 L + m R ) ± y (m^ — m 2 R ) 2 + 4 |X| 


ml = + h 2 t \H°\ 2 , ml = ml + h 2 \H° 2 \ 2 , (33) 

X = \X\ = h t (A t H$ + , (0 < ip < 2vr), (34) 

where mg, ml and At are the soft supersymmetry-breaking mass parameters of the stop sector. 
We assume here fi and A t to be real [], so that (|A|) = (htv 2 /\/Z) \A t + pcotp\ and (e lv ) = 
sign (A t + /rcot/3), but we do not make any assumption on the sign of p and A t . A fourth 
quantity is the mixing angle 9f in the stop sector (0 < Of. < tt/ 2), where: 


sin 2 9i 


2\X\ 

2 2 ~ ’ 

m~ — Tri- 

ii 12 


cos 29^ = 


2 2 
m L — m R 

2 2 ~ ‘ 

ra~ — m~ 

11 C2 


(35) 


Notice that the usual field-independent definition for the angle 9f that diagonalizes the stop mass 
matrix at the minimum, 


sin 2 9i = 


2 m t (A t + pcot/3) 

2 2 

m~ — m~ 

1 1 C2 


cos 261 = 


u 


m Q — m 

2 2 ~ 

771- — m~ 
tl 1 2 


(36) 


leads instead to —7 t/2 < 9{ < ir/2. A fifth quantity can be chosen to be a function of the phase 
difference tp — (p, where ip is the phase in the complex top mass, X = ht H 2 = |X| e ^¥, , with 
0 < p < 27r such that (|X|) = /itu 2 /\/2 and (p>) — 0. A convenient choice is c ^ = cos (ip — <p), 
with (cpp) = ±1. 

The classes of Feynman diagrams that contribute to the two-loop expression of V and affect 
the 0(a 2 ) calculation of the neutral Higgs boson masses [] are shown in Fig. 1. 

Introducing a wave function renormalization (w.f.r.) matrix Z for the Higgs fields, we can 
write 

(AMI) 6 ® = VZ (AAll) 65 ^, (37) 

where 

(AAll) 6 ® = i/qV4 (F; i + 2A II F 3 ) , (38) 

^AAi!^ ^ = h 2 pmt S 29 (F 2 + A^F 2 ) + — h 2 At p s 2e (F% + AA t F% + A^F%) , (39) 

cfF l 

(AA4l) 22 = 2h 2 m 2 Fi + 2h 2 A t mt S 20 (F 2 + AA t F 2 ) +-h, 2 A 2 s 2 q (F% + 2 A^F^) . (40) 

In Eqs. (||)-@), h 2 = 2m 2 t /(v 2 s 2 0 ) is the top Yukawa coupling, and s 2 6» = sin 29j refers to 
the field-independent stop mixing angle defined in Eq. (|36|). The functions F t (i = 1,2,3) are 

5 Notice that Eq. ([d]) implicitly defines our conventions on the sign of the parameter p, and fixes the sign 
in the chargino and neutralino mass matrices. Notice also that our conventions differ by a sign from those of 

6 An explicit expression of the O(aj) corrections to the effective potential can be found in Eq. (D.6) of 
To extend its validity to the general case discussed in the present paper, such equation should be modified in the 
following way: take the couplings A'Lor r and AL. + - - , appearing in the last line of Eq. (D.6); in the corresponding 
dehnitions, Eq. (B.6), multiply the terms m,tS 2 eXt and mt by the factor c v qj. 
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Figure 1: The classes of Feynman diagrams that contribute to the two-loop effective po¬ 

tential and affect the O(af) calculation of the neutral Higgs boson masses [q = (t,b), ip = 


(H,h,G,A,H ± ,G±),h=(h° li2 ,h ± ), q = (i x .t 2 .b L )}. 


written as F{ = Fi + A Fi , where F % are combinations of derivatives of the effective potential with 
respect to the field-dependent quantities, whose explicit expressions are given in Eqs. (28)-(30) 
of 0. The terms A Fi include the renormalization of the common factors multiplying F t (i.e. 
hf, nit, S 20 ), as well as the renormalization of the mass parameters appearing in the one-loop 
parts of the Fj, while A A t Fi and Aencode the renormalization of the factors that are not 
common []. 

We performed the evaluation of (F\ , F 2 , F 3 ) following the strategy of 0 , namely computing 
directly the derivatives of V, without evaluating V explicitly. Before discussing the O(af) two- 
loop terms we recall first, for completeness, the one-loop 0(at) result j 2 ]: 


F^ = 


N c 
16 7 r 2 


■ In 


2 2 

mf mf 

tl t-2 

mf 


F\ l = 


N c 
16 7T 2 


■ In ■ 


m 


m 


f¥ = 


N c 
16 7T 2 


2 - 


77T? + m~ 

-Ji - h. i n . 

2 2 111 

m~ — mf 
1 1 12 


where N c = 3 is a color factor. 

We organize the presentation of our two-loop results as follows. 


(41) 

In this section we give the 


'Indeed, the terms Fi in Eqs. ( fig} ) ( fio| ) are multiplied by different combinations of /r and At. These two 
parameters have different O(at) renormalizations that cannot be absorbed into the Fi, and are then separately 
taken into account. In Ref. |IcJ, which deals with the 0(ata s ) corrections, the expressions analogous to Eqs. (§)- 
© contain neither the w.f.r. Z, nor the terms A ^Fi, since in such case the w.f.r. is not needed and the parameter 
fi does not receive 0(a s ) corrections. In [lo|, what we call here (AAT|) eff and A A t Fi are indicated as (A M%) 
and AFi, respectively. 
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results assuming that all the parameters (masses, VEVs, mixing angles and Yukawa couplings) 
appearing in the one-loop term of Eqs. (|38|)-(41) are expressed in the DR scheme. In the next 
section we give the expressions for the shifts of these parameters that allow to translate the DR 
result into other renormalization schemes. In particular, we specialize the shifts to a convenient 
renormalization scheme, and we comment on alternative choices in the existing literature. 

Since the general analytic expressions of (F \, F 2 , F%) in Eqs. (38) (|f]) are very long, we 
choose not to display them in print. Instead, they can be obtained upon request as a Mathematica 
or Fortran code. Here we present explicit analytical formulae for tuq = mu = Mg, and neglecting 
0(mt/Ms) corrections, but without making any assumption on the size of m A . In this case, the 
stop masses are given, in the gaugeless limit, by rnj^ ^ = Mg + rri 2 ± mt\Xt\, where Xt = 
At + fj.cot/3, and s\ e = 1. Setting consistently the bottom mass to zero, the relevant sbottom 
mass is m~ = Mg. We denote by {F\, F 21 , F'f ) the expressions for the two-loop 0(a 2 ) part of 
the Ft functions obtained assuming that the one-loop part is evaluated in terms of DR quantities. 
In units of h 2 N c /(16 ir 2 ) 2 , we have: 
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where e = (4 —n)/2 and n is the dimension of the space-time. We notice that all the double poles 
1/e 2 generated by two-loop diagrams have cancelled in the functions (i = 1,2,3), whereas 
a few 1/e poles have survived. Also, Q is the DR renormalization scale, Y) = At — fi tan [3, 
Li 2 (z) = — fl dt [ln(l — t)\/t is the dilogarithm function, and 

f 4^~^ Cl 2 (2arcsin yfz) , (0 < 2 < 1) , 

ct>(z) = \ (45) 

{ \ [-4Li 2 (i^) + 21n 2 (^) -ln 2 (42)+7r 2 /3] , (z > 1) , 

where C\ 2 (z) = ImLio(e* 2 ) is the Clausen function and A = yjl — (1 /z). 

Under the same approximations of Eqs. (|42|)-([44|), and in the same units, the 0(a 2 ) contri¬ 
butions to A ^Fi and A A t Fi are: 
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Finally, the w.f.r. matrix Z induces additional contributions to the entries of (A_A/f|) cff . In 
the DR scheme Z is diagonal, ( Z)ij = Zidij , and only Z 2 is different from 1, since we are 
neglecting the gauge couplings and the bottom Yukawa coupling. The entries (1,2) and (2,2) get 
then additional corrections that in units of hf IV 2 /(16 7r 2 ) 2 , and in the same approximations of 


Eqs. (|4^)-(|47D, read: 
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Some comments on Eqs. (12)-' 
the w.f.r. [Eqs. 


are in order. It can be checked that, after the inclusion of 
-(49)], all the leftover divergent contributions to (AAl|) cff do indeed cancel 


out. We have verified that an analogous cancellation takes place also in the general case. The 

matrix (AM|) shows an explicit dependence on the renormalization scale Q, due not only 

to our choice of expressing the various parameters in the DR scheme, but also to the fact that 

the entries of rs(p 2 ) are not physical quantities. Only the solutions to Eq. ([7]) should be Q- 

independent, once the quantities they depend upon are expressed in terms of physical observables. 

Then, to obtain a Q-independent result for the eigenvalues, we should also take into account, in 

2 

the perturbative diagonalization of T s(p 2 ), the contribution induced by (AYf|) :p in Eq. (p4|). 

We have also checked that, for rn A = Ms ^$> m z and neglecting 0(mt/Ms) corrections, 
the O(al) contribution to rnf t due to (AYl|) eff , namely c| (AA4|)^ + 2 sgcp (AY4|)^2 + 
s 2 p (AMD 22 ) gi ven by Eqs. (p^)-([49D, reduces to Eq. (17) of Ref. [ pT| . 


4 Input parameters and renormalization schemes 

The general results of our calculation can be combined with the complete one-loop calculation 
and with the calculation of the 0 (ata s ) two-loop corrections, to give the most accurate determi¬ 
nation of the neutral CP-even Higgs masses available to date. Here and hereafter, when referring 
to the masses of the CP-even Higgs bosons, (m^, m H ), we will always mean their physical (=pole) 
masses. Before presenting our numerical results, we must discuss the input parameters that need 
to be defined beyond the tree level. They are the parameters (m A ,m z , tan (5) that appear already 
in the tree-level result, and the parameters of the top and stop sectors that appear first in the 
O(at) one-loop correction. Other quantities, such as the gluino mass nig, or the strong coupling 
constant a s , appear only at the two-loop level and do not require further specification. The 
same is true for the parameters, such as the SU(2)l x U(1)y gaugino masses (M 2 ,Mi) or the 
remaining sfermion masses, that appear only in the complete one-loop correction but do not 
appear in the 0 (a.ta s + a 2 ) two-loop correction. 
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According to the discussion of Section |], the parameters m A and m z are chosen to be the 
physical masses, while the result for (AAd|) eff presented in Section |3] is expressed in terms of 
DR quantities. We give now the formulae that allow to express our results in terms of input 
parameters given in a different renormalization scheme, which we indicate generically as R. In 
the case of the 0(ata s ) corrections, the procedure was explained in [ 40 ], but we shall complete 
the discussion here to take into full account the 0(ata s + a? ) corrections. 

To obtain the 0(ata s + ck 2 ) correction in the R scheme, we have just to shift the parameters 
appearing in the one-loop term. Indicating, generically, a quantity in the DR scheme as x DR , and 
the same quantity in the R scheme as x, we can write the one-loop relation x = x DR — 5x. Then, 
once the one-loop contribution is evaluated in terms of R quantities, the two-loop 0(ata s + af ) 
corrections in the R scheme can be obtained through: 
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(51) 


(52) 


A At F t = A At Fi + 6 -^ F t u , A^Fi = , (i = 2,3), (53) 

At /I 

where all the quantities that appear in Eqs. (|50|)—(|53|) are meant in the R scheme. We notice 
that the shifts appearing in Eqs. (|5^) ( |53|) are not all independent, because of the relation (|36|). 

We choose to express our result identifying mt with the pole top mass, while the top Yukawa 
coupling will be defined indirectly by ht = v^2 mt/{ysp). For the electroweak symmetry-breaking 
parameter v, we will use the value obtained in terms of the precisely known muon decay constant 
Gfj, according to the relation v = (y/2 G ^) -1 / 2 = 246.218 GeV. These choices set Srrit and 5v to 
be 


5m t = R et t (m t ), 


5v_ _ 0 ) 

v 2 m^r 


(54) 


where St(mt) is the finite part of the top quark self-energy, evaluated at an external momentum 
equal to the top mass, and is the finite transverse part of the W-boson self-energy, 

evaluated at zero external momentum. 

The parameter tan /3 appears in the tree-level result, in the stop mass matrix and also, 
implicitly, in ht. We recall that, at the tree level and in the gaugeless limit, the stop masses 
and mixing angle depend on mt and five more parameters, (mg, mu, At, /t, tan (3). The last three 
enter the stop mass matrix only in the combination Xt, so that (m^ , , s 2 q) are completely 

defined in terms of the three quantities (mQ,mjj,Xt). Yet, the dependence of (AA4|) cfr on 
(At, p, tan/3) is not through the combination Xt but is more complicated, as can be seen from 
In the case of the 0(at.a s ) corrections, since p and tan/3 do not receive G(a s ) 


Eqs. (HM 

one-loop corrections, the specification of the three parameters of the stop sector, together with 
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nit, is actually sufficient to fully determine (A.Mg) cfr . The situation is different for the O(af) 
corrections. Since A t , /z and tan/3 are all corrected by 0(at) one-loop effects, their individual 
definition is required, in addition to the specification of (m;, v. rri^. rrif 2 ). However, because of 
the relation (36), we can actually trade one of the parameters (At,/ a, tan/3) for s 2 e»- 

While there is a well known OS definition for the masses, and an OS definition for the stop 
mixing angle can be also conceived, it is not clear what meaning should be assigned to an OS 
definition of (A t , /qtan/3). For instance, they could be related to specific physical amplitudes. 
However, given our present ignorance of any supersymmetric effect, such a choice does not 
seem particularly useful. Then, it may seem simpler to assign all the parameters related with 
the stop sector, i.e. (rriQ, mjj, A t , /z, tan/3), in the DR scheme, at a chosen reference scale Qq. 
However, as noticed in flioll , in the DR scheme the limit of heavy gluino is very violent, since 
the 0(ata s ) corrections to the Higgs masses include terms proportional to rrig and m~. This 
powerlike behavior is actually not present if the stop masses and mixing angle are assigned in a 
suitable OS scheme, since it is cancelled by the relevant shifts in Eqs. (|50|)— (53). In fact, with the 
latter choice the 0(ata s ) corrections grow only logarithmically as rrig increases. In this situation, 
we think that the use of OS stop masses and mixing angle should be preferred, while for the 
quantities /z and tan (3 we are still going to employ the DR definition, at a reference scale that we 
choose near the present central value of the top quark mass, Qq = 175 GeV. In this framework, 
we treat At as a derived quantity, obtained through (|36|), whose shift is given by 
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-<5m| | 5s 2e 5m t 

S20 rnt 


2 2 
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11 12 
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(55) 


According to the above discussion, we are going to take <5/z = 5/3 = 0. The requirement of 


OS stop masses specifies Sm~ and 6m~^ uniquely as 


bm\ x = Re fin (m|), 


5m h = ReII 22 (™| 2 ), 


(56) 


where iijHm- ) (i = 1,2) are the finite parts of the diagonal stop self-energies, in the mass 
eigenstate basis, evaluated at external momenta equal to the corresponding masses. For the stop 
mixing angle, several OS prescriptions are present in the literature (for a discussion see, e.g., 


[12, [l3], 14] and references therein). We find the following ‘symmetrical’ definition (lj, H3| 
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(57) 


such that 5s2o/s2e = 2 cot 26f 56 suitable for our O(af) calculation, while the choice 


56- t 


n 12 (m|) 

2 2 ’ 
m~ — m~ 

11 C2 


(58) 


employed in || for the calculation of the 0(ata s ) corrections, seems less justified at O(<%%), since 
the divergent parts in the ‘bare’ version of Eq. (|5^) do not match. In Eqs. (|57|)-(^1), fli 2 (p 2 ) 
stands for the off-diagonal stop self-energy. 

Explicit formulae for the shifts relevant to our calculation are provided in the Appendix. 
These formulae include only the contributions that come from the top-bottom sector of the 
MSSM, and are specialized to the gaugeless limit we are considering in this work. General 
expressions for the top and stop self-energies at one loop can be found in the literature |5], [h| . 
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5 Numerical results 


In this section we discuss the effect of our O(of) two-loop corrections on the masses of the 
CP-even Higgs bosons, nih and rn H . As anticipated in the previous section, we take as input 
quantities the top pole mass m* = 175 GeV, the electroweak symmetry-breaking parameter 
v = 246.218 GeV, the OS masses (m A , m^, m^ 2 ), the OS stop mixing angle 6j. and the DR 
quantities (fj,, tan /?), evaluated at the reference scale Q o = 175 GeV. To facilitate the comparison 
with the existing analyses of the two-loop corrected Higgs masses, we also refer to the unphysical 
parameters (mg S ,m° s ,T ( os ) that can be derived by rotating the diagonal matrix of the OS stop 
masses by an angle However, it must be kept in mind that such parameters are ‘on-shell’ 
only in the sense that they are extracted from the OS masses and mixing angle, but they do not 
correspond directly to any physical quantity. The gluino mass rrig and the left sbottom mass 
rrtr enter only in the two-loop part of the calculation, thus their precise definition amounts to 
to a three-loop effect and is irrelevant to the perturbative order we are working at. Anyway, 
is not really a free parameter, i.e. it has to be derived from the parameters of the stop sector: 
a simple choice consistent with our approximations is = rriQ S . Finally, the remaining input 
quantities are m z = 91.187 GeV and the strong coupling constant a s (m z ) = 0.118. 

Since we are interested here in pointing out the relevance of the O(a^) corrections with 
respect to those 0(at) and 0(ata s ), we include in the one-loop part of the calculation only the 
contribution of the top-stop sector of the MSSM, as given in ||. A general discussion of the 
precise determination of the Higgs masses as functions of the whole set of MSSM parameters 
will be presented elsewhere. For the two-loop 0(ata s ) corrections we use the simple analytical 
formulae of [|H]] , equivalent to the results of the diagrammatic calculation |8| implemented in the 
Fortran code FeynHiggs |lG[], with the advantage of being simpler and more flexible for the choice 
of the renormalization scheme. For the O[o^) corrections to the CP-even Higgs mass matrix, 
we use our complete analytical formulae in the version appropriate for OS input parameters. 
In the analysis of [|j], some ‘leading logarithmic’ O(a^) corrections have been added to (A4|) 22 , 
proportional to powers of In (m~ mf /mf) and obtained by renormalization group methods fl. In 
the following we will compare our complete O(af) result with the renormalization group result, 
and we will point out that, for some choices of the SUSY parameters, the logarithmic corrections 
amount only to a fraction of the full ones. 

Fig. H shows TYiji as a function of m A , for tan/? = 2 or 20 and A ® 3 = 0 or 2 TeV. The 
other input parameters are chosen as m,Q S = m ® 3 = 1 TeV, n = 200 GeV, nig = 800 GeV. The 
cases with A t os = 0 and A f ° 3 = 2 TeV correspond, respectively, to the so-called ‘no-mixing’ 
and l mi ,-max’ benchmark scenarios considered in the experimental analyses |Gj. The values 
of the stop masses and mixing angle are approximately, neglecting small D-term contributions, 
mi = mi 2 = 1015 GeV, = 0 in the no-mixing case, and = 1175 GeV, m ^ 2 = 825 GeV, 
9i = | in the case X® s = 2 TeV. The curves in Fig. |2] (and also in the other figures, unless 
indicated differently) correspond to the two-loop corrected Higgs mass at 0(ata s ) (short-dashed 
line), at 0 (cqa s + af) including only the logarithmic corrections of [|J (dot-dashed line), and at 
O(ata s +o.t ) including our full computation (solid line). The one-loop result for m/j is also shown 
for comparison (long-dashed line). It can be seen from the figure that the 0(a,ta s ) corrections 
are in general large, reducing the one-loop result for mh by 10-20 GeV. On the other hand, the 
O(ctt) corrections tend to increase mh■ for small stop mixing they are generally small (less than 
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X os = 0 


X os = 2 TeV 




m A (GeV) m A (GeV) 

Figure 2: The mass rrih as a function of m A , for tan (3 = 2 or 20 and X^ s = 0 or 2 TeV. The 
other parameters are m,Q S = rrijj S = 1 TeV, fi = 200 GeV, m,g = 800 GeV. The meaning of the 
different curves is explained in the text. 


tan[3 = 2 tanP = 20 



Figure 3: The masses ( mh,m H ) as functions of m A , for tan/3 = 2 or 20 and V f os = 2 TeV. The 
other input parameters are as in Fig. |j. 
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2-3 GeV), whereas for large stop mixing they can reach 7-8 GeV, i.e. a non-negligible fraction 
of the 0(ata s ) ones. Moreover, it appears that the ‘non-leading’ corrections included in our 
full result are always comparable in size with the ‘leading logarithmic’ ones, and are even more 
significant than the latter in the case of large stop mixing. For m A 3> m z , we find good numerical 
agreement with the value of rrih obtained from the approximate formula (20) of Ref. jll| |. 

In Fig. H we plot the one-loop and two-loop corrected values of both nih and m H , in the 
region m A < 200 GeV, for the case of large stop mixing. The other parameters are chosen as 
in Fig. |2|. It can be seen that the radiative corrections affect significantly m H only for quite low 
values of m A (m A ^ 200 GeV for small tan/3 and rn A < 140 GeV for large tan/3). We stress that 
the determination of both mass eigenstates is made possible by our implementation of the EPA, 
whereas other existing two-loop calculations in the EPA i 10 addressed only the corrections 
to mh. in the limit m A rn z . 

Fig. U shows the corrections to the effective CP-even Higgs mixing angle a, for the case 
Xf s = 2 TeV, in the region of moderately low m A (for large m A the mixing angle sticks to 
its tree-level value, /3 — ^). The other input parameters are as in Fig. |2|. We have chosen to 
plot the combination sin 2 (/3 — a) that enters the cross-section for the process e + e~ —► hZ. 
Again, it turns out that the 0(a 2 ) corrections can reach 40% of the 0(ata s ) ones, and that the 
renormalization group method provides only a fraction of the full two-loop Yukawa corrections. 
In the interpretation of Fig. [|, however, the reader should keep in mind that our renormalized 
angle a does not exhaust the 0 (ata s + a 2 ) corrections to the considered process. 

In Fig. [| we explore in more detail the dependence of the two-loop corrected Higgs masses 
on the stop mixing parameter, X® s . We set tan/3 = 2 and m A = 120 GeV or 1 TeV. Taking 
advantage of the fact that our analytical formulae are valid for arbitrary values of the MSSM 
parameters, we split the values of the diagonal entries of the stop mass matrix, choosing m.Q S = 1 
TeV and m^ s = 700 GeV (in fact the choice mg = mu , commonly considered for simplicity, 
seems to us quite unnatural, since it is not preserved by potentially large O(at) radiative correc¬ 
tions). Once more, it can be seen from Fig. |5] that, for large stop mixing, the 0(a 2 ) corrections 
to mh (and also to m H , for small m A ) amount to a significant fraction of the O(ota s ) ones, and 
are not well accounted for by the ‘leading logarithmic’ terms. 

Finally, in Fig. we show the variation of our result for r% under a change of renormalization 
scheme. We use the same input parameters as in Fig. ||, for the case m A = 1 TeV. The solid line 
in Fig. d corresponds to evaluated with our choice of OS input parameters, as described in 
Section 4. For comparison, we have converted the OS input parameters into DR quantities by 
means of Eqs. (|54|)— (57), and we have used the latter as input of the pure DR formulae, using 
the original F 2( instead of the shifted Ff e of Eqs. ([5(]) (|53|) . The results for rrih are plotted as 
a dashed line. The differences between the two curves can be interpreted in part as three-loop 
effects, induced by the variation of the parameters that enter in the two-loop formulae, and 
in part as two-loop 0(g 2 at) effects, due to the fact that the counterterms of Eqs. (|50| ) ( |53|) 
account only for the scheme dependence of the 0(at ) EPA corrections, whereas in the one-loop 
corrections we include also the 0(g 2 ) D-term contributions and the momentum dependent parts 
of the propagators. The figure shows that for moderate V f os there is good agreement between 
the two curves, the differences being 0(1 GeV) or smaller. On the other hand, for Y t os large 
and negative the DR calculation of mh blows up. The reason is that, for such values of Y f os , 
the stop masses and mixing angle are subject to very large O(at) corrections proportional to A 2 . 
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tan(3 = 2 tan(3 = 20 



m A (GeV) m A (GeV) 

Figure 4: The effective CP-even Higgs mixing angle a, in the combination sin 2 (/3 — a), as a 
function of m, A , for tan/3 = 2 or 20 and ?Bq S = m® s = 1 TeV, X t os = 2 TeV, /i = 200 GeV, 
m,g = 800 GeV. The meaning of the curves, explained in the text, is the same in the two frames. 



X os (TeV) 


X os (TeV) 


Figure 5: The CP-even Higgs masses as functions of the stop mixing parameter X ^ s , for m A = 
120 GeV or 1 TeV. For large m A only rn^ is shown. The other input parameters are = 1 
TeV, = 700 GeV, tan (3 = 2, fi = 200 GeV, nig = 800 GeV. 
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X OS (TeV) 


Figure 6: The mass as a function of the stop mixing parameter . The input parameters 
are chosen as in the right frame of Fig. |5[ The solid line shows m/j evaluated with OS parameters, 
whereas the dashed line shows m/,. evaluated with the corresponding DR parameters. 

As discussed in the previous section, a similar phenomenon would also occur for a large gluino 
mass, because of the non-decoupling properties of the DR renormalization scheme. 

6 Conclusions 

In this paper we computed the 0(a 2 ) two-loop corrections to the mass matrix of the neutral CP- 
even Higgs bosons in the MSSM, extending earlier results 10] on rrih. Using the formalism of 
the effective potential, we obtained complete analytic expressions for the momentum-independent 
part of these corrections, valid for arbitrary values of rn A and of the parameters in the stop sector. 
The inclusion of momentum-dependent effects was explained in Section 2. Our results can be 
efficiently inserted in general computer codes, to obtain 0(ata s + a 2 ) corrected expressions for 
the masses ( rrih,m H ). 

Our general analytic expressions are too long to be useful if explicitly written on paper, thus 
we decided to make them available as a Fortran or Mathematica package. As an illustration, 
we wrote down explicit analytical formulae for the simplified case where tuq = mu = Ms and 
0(mt/Ms ) corrections are neglected, but m A and the stop mixing parameters are left arbitrary. 

We have presented our results in such a way that the input parameters of the top and stop 
sectors can be assigned either in the DR scheme or in any other renormalization scheme, in 
particular those with on-shell definitions of the masses and of the stop mixing angle. We have 
also included in our computer package the analytic result for the 0(a 2 ) two-loop corrected 
relation between m 2 , defined in the DR scheme, and the EPA mass m 2 A , with the remaining 
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parameters assigned in the DR scheme. 

As discussed in Section ||, the 0(a 2 ) corrections have a significant impact on the predictions 
for rrih and m H . While the 0(ata s ) corrections typically decrease the values of (m^, m#), the 
0 (a 2 ) ones typically increase them, producing a partial compensation: in the case of large stop 
mixing, the latter can reach up to 40 % of the former, and shift mh upwards by 7-8 GeV. Also, 
our general result for the O(af) corrections, both to the masses ( mh,m H ) and to the mixing 
angle a, is quite different from the one obtained via a ‘leading logarithm’ renormalization group 
analysis: the shift in mh is twice as large for small stop mixing, about three times larger for large 
stop mixing. 

The present paper and its companion, Ref. [10], contain the most advanced calculation (so 
far) of the two-loop corrections to the MSSM neutral Higgs boson masses. The 0(ata s + a 2 ) 
corrections to m 2 H are still affected by some residual uncertainties for large m A , because of 
momentum-dependent effects not computed so far. However, the relative importance of such 
corrections is small, since the tree-level value of m H grows with m A . As explained in Section 2, no 
such uncertainties exist for mh- The main residual two-loop corrections to the Higgs masses are 
controlled by the electroweak gauge couplings and by the bottom Yukawa coupling. In particular, 
two-loop corrections proportional to atb are expected to be numerically non-negligible only for 
tan/? > 40, i.e. when is comparable with at, despite the fact that mi, <C m t . The 0(aba s m 2 ) 
and 0 {a 2 m 2 ) corrections could be computed by performing simple substitutions in the formulae 
for the top case, whereas the 0(ata},m 2 ) and 0(atabm 2 ) corrections require further work. We 
also recall that the bottom mass mj can receive quite large threshold corrections [18] that could 
induce additional complications. 

In conclusion, our work should lead to a more accurate interpretation of the experimental 
searches for the neutral MSSM Higgs bosons at LEP, the Tevatron, the LHC and other possible 
future colliders. The importance of the new O(o'f ) effects we have computed will increase fur¬ 
ther when the top quark mass is measured more precisely (we recall that the present CDF/DO 
average [ 191 is mt = 174.3 ±5.1 GeV). Then, we can hope for the next step, the discovery of 
supersymmetric particles and supersymmetric Higgs bosons. 
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Appendix: One-loop self-energies 

We provide here the explicit formulae for the one-loop self-energies of the top quark and squarks 
and of the W-boson. As discussed in Section 4, these formulae are required to relate the pa¬ 
rameters mt, rnj , mj , S 20 and v defined in the DR scheme to the corresponding quantities 
chosen in that Section. Since we are only interested in the terms that generate 0(a 2 ) corrections 
at two-loops, we neglect all gauge couplings and we keep ht as the only non-vanishing Yukawa 
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coupling. This leads to the simplified mass spectrum described at the beginning of Section 3. 
The real finite parts of the self-energies for the top and the stops are: 


Y t (m t ) = 


h 2 t „ j (-, 5 2 V , 3c /3 f, , rn 2 A \ 3 p 2 / 1 , /i 2 

32 7 T 2 mt | V 2 S Pj ln Q 2+ 2 rnj [ n Q 2 ) 2 m 2 ( n Q 2 


4 4 ™ 2 


+5 5^ - 1 + — I 1 - j B o (m t ,0,m A ) + c p I 2 - —f I B 0 (rn t , rn t , rn A ) 


1 T rri~ ( rnj \ rnj - mj + p 2 , 

+ 2 4 1 - lD # + 4 ‘ 


+ (4 4) + (4 ~> 4) 


nil (m?) = { 4 ( x + 4) ^0 (mi) - A 0 (m|) - (1 + s 2 0 ) A 0 {p 2 ) 

+4 n 0 (m~ L ) + (cj g - Nc 2 1 A 0 (m|) + ^£±1 S \ Q a 0 ( m ? ) 

+(mf x - m 2 - 4) ^(m? , m|, /r 2 ) + s| (m| - /r 2 )S 0 (m| i , 0, /x 2 ) 

s 2 p (2 m t + s 20 At) 2 B 0 {m 2 x , m|, 0) + c| (2 m t + s 20 T t ) 2 B^m 2 ^, m|, m A ) 

+ s /3 (1 + c 2 e) A t 2 ^o( m f 1 , m| 2 , 0 ) + c| (1 + c|e) T t 2 Bojm ^, m| 2 , mi) 

+ s 2 p (m t eg + X t Sg) 2 Bq (m^, m^, 0) + 4 (m t eg + Y t s e ) 2 Boim 2 ^, mj^, mi) | , 

(A2) 


ni 2 (p 2 ) = ^2 { s 20 (p 2 - p 2 )B 0 (p 2 ,0,p 2 ) - s 29 A 0 {p 2 ) 

+s 2 e 4 ^0 (m 2 ) + s 20 A 0 (m^) + (1V C + 1) c 26 s 28 A 0 (m|) - A 0 (m|) 

+s| c 2e A t (2 m t + s 2 g X t ) 4 0 (p 2 , m|, 0) + c| c 2e 1) (2 m t + s 2e 4) 4 0 (p 2 , m|, mi) 

+s|c 2e A t (2m t - s 20 X t )S o (p 2 ,m|,O) +c|c 2e l) (2m t - s 2 gY t ) B 0 (p 2 ,m|,mi) 
-2 (mt ce> + ^ s 0 ) (mt se - A t eg) B 0 (p 2 , m-^, 0) 

-24 (mt c 0 + 4 «e) (m t s 0 - 4 Q;) 4 0 (p 2 , m^, m 2 ) j , (A3) 
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where 


Ao(m 2 ) = m 2 1 — In 


m " 


and -Bo denotes the real finite part of the Bq Passarino-Veltman function, i.e. : 

Jo Q 


(A4) 


(A5) 


An explicit expression for Bq can be found e.g. in [20]. The function II 22 (m~ ) can be obtained 

4_1J Z2 ' 


from the right hand side of Eq. ( |A2| ) with the replacements m~ m~ , eg —> sg , sg —>■ —eg 
(also implying s 2 e ->■ -s 2 0 ,c 2( ? -> ~c 2 g). 

Finally the identification v = (\/2 GL) -1 / 2 = 246.218 GeV requires the self-energy of the IT 
boson at zero momentum transfer. Including only the contributions from the top-bottom sector 
relevant for our calculation: 
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